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Abstract 

The branching ratios for the neutral and charged B decay channels to rf K and r\K* are 
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both these discrepancies, without jeopardizing those results which are in agreement with the 
Standard Model. We found that it is possible to have a parameter space satisfying all these 
requirements, including that of a low strong phase difference compatible with the color trans- 
parency argument. Furthermore, we put a robust bound on the relevant coupling, which is two 
orders of magnitude better than the existing one. 
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1 Introduction 



The flavour structure of the Standard Model (SM) is, arguably, the best proof that there is new 
physics lurking above the electroweak scale. We do not know how this new physics should manifest 
itself, so one is forced to consider all possible ways to uncover it. Moreover, any experimental result 
that deviates from the SM expectation is bound to receive close scrutiny for new physics. Apart 
from the direct production of new particles, one can look for their indirect effects, i.e., how the 
operators at a high energy affect the low-energy observables. In this respect, the e + e~ B factories 
at Cornell, SLAC, and KEK are doing a commendable job by churning out a huge amount of data 
on various B decay modes, including the branching ratios (BR) and CP asymmetries. 

Even at this stage, there exist several hints that everything is not as one should expect from 
the SM, if we take the experiments seriously. The SM is in no way ruled out, since there is much 
scope for theoretical uncertainty in low-energy QCD, but it will be under considerable strain if 
the experimental results persist over the next few years. In a large set of data, most of which 
are in perfect agreement with the SM, there are three so-called sore thumbs sticking out: (i) The 
direct and mixing-induced CP-asymmetries for the mode B — > 7r + 7r~ where Belle and BaBar are 
not consistent [1, 2], (ii) The abnormally high branching ratios for the B — > rj'K and B — > r]K* 
modes, incompatible with the SM prediction with factorization [3], and (iii) The measured value 
of sin 2/3 from B — > 4>Kg [4, 5] 1 . (Last two points are closely related; see, e.g., [7].) Let us briefly 
comment on these anomalies. But first we wish to make our notation clear: we use B as a general 
shorthand for both B° and B°, and also for the charged B mesons where the mention of charge is 
either unnecessary or self-evident, and the BRs are assumed to be averaged over the CP conjugate 
states. Moreover, we use the symbol r\K to denote the final states with either an rj or rf and a K 
or K* (neutral or charged). 

The first anomaly exists only for Belle — BaBar data is completely consistent with the SM 
expectation, if one uses perturbative QCD calculations [8] or QCD-improved factorization model 
[9]. However, if we take the Belle data seriously, existence of new physics is hinted. A possible 
solution is discussed in [10]. 

The second and third anomalies are the main theme of this paper. Before one confronts the 
data with any specific model, three things are obvious. First, even if there is new physics, it is very 
much flavour-sensitive, since modes with the same topologies are not equally affected. Secondly, 
the effect of new physics is of the order of unity, considering the change in the BRs of B — > r\K 
modes, and the CP asymmetry from the B — > <f>Ks mode. Thus, the effective operators generated 
by the new physics should better be tree-level at the high scale, unless one prefers a strong-coupling 
theory. Thirdly, the new physics amplitude should have weak (and possibly strong) phases different 
from the corresponding SM amplitude to generate the observable CP asymmetries, particularly in 
the B -> 4>K S channel [11]. 

These constraints severely limit the choice of plausible new physics options 2 . Among those 
that survive, one of the most respectable options is R-parity violating supersymmetry (RPV) 
[13, 14]. RPV automatically generates flavour-specific tree- level operators which can leave their 
signals in some, but not all, B decay modes. 

1 There are some other anomalies, e.g., lack of a consistent fit of the BRs of the B — » PV channel [6]. It is not 
clear at this moment whether this is a feature of QCD factorization models, which was used in the analysis, or is 
true in general. 

2 There may be new physics with loop-induced and/or flavour-blind operators, but their effects will be very hard 
to detect in the B factories. However, some supersymmetric resolutions of the <j>K puzzle have been proposed [12]. 
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The effect of RPV on B decays has been extensively discussed in the literature [10, 15, 16]. 
It has been proposed as a solution to the B — > <j)Ks CP asymmetry puzzle [17], and even a 
simultaneous solution to both the rjK and <j)Ks anomalies has been proposed [18]. 

What is new, then, in this paper? There are two main points where our analysis differs 
from [18]. This is the first time that one takes all the relevant parameters, including the final- 
state strong phases, into account, and tries to explain all available data, including the direct CP 
asymmetry results. Secondly, we have not looked at isolated points in the parameter space as 
possible solutions, but have made a complete scan over the whole parameter space, and find out 
the regions which satisfy all the existing data. In fact, we introduce only one new physics operator, 
which, we find, does the fitting quite well. Moreover, we get a stronger bound on the relevant 
coupling, two orders of magnitude better than the existing one. 

The discrepancies from the SM expectations show up in the BRs of different B — > r/K modes 
(see our convention before), and the value of sin(2/3) extracted from B — ► <pKs decay. The 
numerical values are given in Section 2. By itself, the CP asymmetry data has nothing wrong in 
it, but it is definitely not consistent with the value of sin(2/3) extracted from B — > J/ipKs (and 
other charmonium channels) even at 2a. 

That the experimental numbers are not consistent with the SM predictions for rfK and r/K* 
modes are known for a long time. Many solutions for this are proposed, ranging from introducing 
a charm content of rj' or anomalous rj' coupling to gluons to the introduction of new physics like 
RPV [15, 19]. For our calculation we take the charm contents in both these mesons to be zero. To 
have an estimate of the SM BRs for these modes, we assume naive factorization (NF) to hold good. 
Note that the BRs are more or less stable with the variation of the number of effective colors, N c , 
which takes into account the nonfactorizable effects [20]. This, however, is not true for the generic 
B — > <j)K modes, which is a penguin and the BR depends sensitively on N c , falling quickly as N c 
increases. Thus, if one takes a fixed value of N c , the BR predictions (within the context of the 
SM) for the r)K* and rfK modes are sort of reliable, but that is not the case for the <j)K mode. 
In fact, it has been shown that dynamical enhancement of the penguin operators may lead to a 
twofold increase of the BR over the expectation from the NF model [21]. Keeping in mind such 
sensitivity which results from a numerical cancellation between different Wilson coefficients (WC), 
we think it prudent to present the analysis both by including Br(B — ► (pK) as a constraint and 
by discarding it. We will show that the tighter constraints on the parameter space come from the 
CP asymmetry data and the B — > rjK BRs, and relaxation of the B — > <fiK BR constraint does 
not affect the solutions in any perceptible way. 

The paper is arranged as follows. In the next section we enlist all our input parameters, rele- 
vant formulae and experimental data. In Section 3, we introduce the RPV effective hamiltonian, 
and show how it affects the modes we are interested in. Section 4 deals with the numerical results, 
and we display the allowed regions of the parameter space that satisfy all data. In Section 5 we 
summarize and conclude. 
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2 Input Parameters 



We consider only those modes which are, even in the SM, governed by the b — > sqq transitions, 
with q = u, d or s. The relevant effective four-Fermi hamiltonian reads 



K 



eff 



Gf 
V2 



10 



VvbKs E c i (»)o i (v l )-v fb v£ y E i c i (ji,)o i (ri 



i=l,2 



i=3 



(1) 



where the operators Oj has the standard form 

Ox = (s a ^P L up)(up^P L b a ), 

03,5 = {Sa^Plba) E fa$Tfl* P L,RQl3), 
q=u,d,s 

3 

07,9 = ^(s a ^P L b a ) {e q qp^ P R,LQl3)i 

q=u,d,s 



(2) 



with a and (3 being the color indices, and e q the charge of the corresponding quark. The opera- 
tors Om are obtained from Om-\ by color-singlet <-> color-octet transformation. The projection 
operators are Pl(Pr) = 1 — (+)75- 

The effective WCs for the transition 6 — > s are evaluated at the scale ji = mft/2 at next-to- 
leading-log (NLL) precision in naive dimensional regularization (NDR) scheme, with m t = 170 
GeV, a s {m z ) = 0.118, a{m z ) = 1/128, and the QCD scale A^L = 225 MeV. They typically 
change by 10% if we vary the QCD scale by about 60 MeV. The values, which are taken from [20], 
read 



d = -0.33, C 2 = 1.16, C 3 = 0.022 + 0.003i, C 4 = -0.051 - 0.009i, 
C 5 = 0.016 + 0.003i, C 6 = -0.063 - 0.009i, C 7 = -(1.2 + 1.3s) x 10~ 4 , 
C 8 = 5 x 10" 4 , C 9 = -(101 + 1.3i) x 10~ 4 , C 10 = 20 x 10~ 4 . 



(3) 



At \x = rrib/2, the current quark masses (in GeV) are taken to be 

m u = 0.0042, m d = 0.0076, m s = 0.122, m c = 1.5, m b = 4.88. (4) 

The masses for the mesons B°, B~ , tt, rj, rf, K and K* are the corresponding central values as 
given in [3]. 

The meson decay constants (in GeV) are: 

U = 0.133, f K = 0.158, f K , = 0.214, U = 0.233. (5) 
The rj and rj' decay constants are obtained from f m 



fU = fvs cos 6 fm sin 



_ f m sm9 f m cos 

J if — 7= 1 



y/6 



f s = - 

Jij 
frf = ~~ 



0.214, /^ = 0.233. 
1.10/* and f m = 1MU [15]: 
2f m cos f m sin 6> 



V6 



sing /^ cos 



\/6 



V3 ' 



(6) 



The mixing angle is taken to be —22°. Thus the numerical values are = 0.099, /* = —0.103, 
f™, = 0.051, /*, = 0.133. This shows why the strange quark plays a dominant role in decays 
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involving an rj '. These values differ slightly from those given in [20] using a two- angle mixing 
scheme taking into account the coupling of gluons to rj and rf; the corresponding numbers are 
0.077, —0.112, 0.063 and 0.141 respectively. Our results do not change appreciably if we use the 
latter set. 

The magnitude of the CKM elements are taken from [22] which uses a fit based only on the 
unitarity of the mixing matrix. The error limits are at 95% CL. 



We could even have used the numbers from the standard CKM fits with Amg as one of the inputs, 
since, as we will show later, the RPV couplings considered here do not affect the B° — B° mixing 
amplitude. This, however, is not true in general; e.g., see ref. [10]. 

The transition formfactors [23] at q 2 = are given by 



F (B -► K) = 0.38; F (B -► rj) = 0.145; F (B rj') = 0.135; A (B K*) = 0.32, (8) 



and -Fo(0) = F\(0). One could have used the so-called 'hybrid' formfactors using both lattice 
QCD and light-cone QCD [20]. However, they are completely consistent with the numbers quoted 
above, and the B — > rjK BRs are only mildly affected by the latter choice. Moreover, the SM 
spread in these BRs is more or less taken care of by a scan over N c . 

The CP asymmetry for B — ► J/tpKs is not modified by our choice of RPV parameters (more 
on this in the next section). This helps us to take the SM value of the angle f3 to be given by [24] 



We show all our results with the central value of sin(2/3), since the uncertainty has negligible effect 
over the final results. 

For the rjK modes, the data reads [3, 6, 25, 26]: 



\V ud \ 

\Vcd\ 

\Vtd\ 




(7) 



sin(2/3) = 0.79 ±0.10. 



(9) 



Br(B + rj'K + 
Br(B + -► r]K* + 
Br(B° -► r]'K° 
Br(B° 7]K*° 
Br(B + r]K + 
Br{B + 7]'K* + 
Br(B° r]K° 
Br(B° r]'K*° 
A CP (B ± -» rj'K* 
A CP (B ± rlK^ 
Af P (B -+ n'K s 
A 1 Q i p(B -+ n'Ks 



= (75 ± 7) x 10" 6 

= (25.4 ± 5.6) x 10" 6 

= {58±\j) x 10" 6 

= (16.41 ±3.21) x 10 

< 6.9 x 10~ 6 



,-6 



< 35 x 10 



,-6 



< 9.3 x 10 



,-6 



< 24 x 10 



,-6 



= 0.11 ± 0.11 ± 0.02 (BaBar) 
= 0.015 ± 0.070 ± 0.009 (Belle) 




(10) 
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where 



A CP (B 



1-|A| 2 
1 + |AP' 





(11) 



with 



A 



(12) 



being the mixing phase in the B° — B° box (2/3 in the SM). Note that BaBar uses a convention 
which differs by a minus sign from our convention of Acp (see eq. (1) of [25]), and the convention 
of Belle differs by a minus sign too, in both direct and mixing-induced CP asymmetries. Also note 
that if B — > rj'Ks were a pure 6 — > s penguin, the mixing-induced CP asymmetry would have been 
sin(2/3); but the presence of a tree- level 6 — ► uus term makes the calculation more complicated. 
That is why we have used the data from [26] which gives directly the direct and mixing-induced 
CP asymmetries, rather than [5], which quotes an effective value of sin(2/3) = 0.76 ± 0.36. 

Among these data, the measured BRs are definitely higher than the theoretical predictions; 
the discrepancy is about a factor of 2 to 3 for the rj'K modes and almost a factor of 10 for the 
r/K* modes within the NF model (this depends on the choice of the regularization scale, and the 
effective number of colors). The expressions for the amplitudes can be found in the appendix, 
and the numerical values of the theoretically expected BRs in tables 8 and 10, of [20]. However, 
there are other models, as the perturbative QCD [21], which predict a dynamical enhancement 
of the penguin amplitude due to the higher- twist corrections. Judging by the enhancement one 
may get in these models for charmless modes, one expects at most to gain a factor of 2 over the 
factorization amplitude. This, evidently, does not serve our purpose. Let us just point out again 
that we take neither any charm content in 77 or 7/ nor any anomalous rj(rj') coupling with gluon. 
The CP asymmetry data are consistent with the SM prediction; if one assumes these decays to be 
penguin dominated as a first approximation, the direct CP asymmetry should be zero, whereas the 
mixing-induced CP asymmetry should be just sin(2/3). The tree-level charged current operators 
may change that prediction. 

The modes B — ► 4>K (both neutral and charged) have also been measured, and the averaged 



where the near equality is expected from the isospin symmetry. The direct CP asymmetry is 
measured by BaBar [25] for the charged mode and by Belle [5] for the neutral mode: 



Note again that both BaBar and Belle conventions differ from ours in an extra minus sign. The 
mixing-induced CP asymmetry, which should give sin(2/3) to a very good approximation, is (at 



BRs are [6]: 



Br{B + -► (j)K + ) = (8.58 ± 1.24) x 10~ 6 , 
Br(B° -► <f)K°) = (8.72 ± 1.37) x 10~ 6 , 



(13) 



A CP {B ± ^cpK ± ) = 0.05 ±0.20 ±0.03 
A CP {B 4>K S ) = 0.56 ±0.41 ±0.12. 



(14) 



la) [4, 5]: 



sm(2p) B ^ Ks = -0.19±g;|g ± 0.09 (BaBar) 
sm(2(3) B ^ Ks = -0.73 ±0.64 ±0.18 (Belle). 



(15) 
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Though their central values are not compatible, one may still make an average: 

am(20) ave = -0.39 ± 0.41. (16) 

We use the averaged values of BR and mixing-induced CP asymmetries as our input parameters. 
For the direct CP asymmetry, we use the Belle data, but it may be observed later that the in the 
entire allowed parameter space the direct CP asymmetry is rather small (between 0.13 and 0.25) 
so that the BaBar numbers are not in trouble. 

The BR for the B — > (f>K mode is a sensitive function of N c . It has also been shown in [21] 
that one needs to take into account the annihilation and nonfactorizable contributions, which push 
the BR up by almost a factor of 2. While the BR is in perfect agreement with the NF model, 
there is no model which can explain the data on mixing-induced CP asymmetry without invoking 
new physics. For this purpose, in our analysis, we use the NF model to calculate the BR keeping 
N c a free parameter, which we take to be the same for both cf>K and r\K modes just for simplicity. 
We will also comment on what happens when one relaxes the BR constraint on <pK. 



3 R-parity Violating Supersymmetry 

R-parity is a global quantum number, defined as (_i) 3B +- L + 25 ' j which is +1 for all particles and —1 
for all superparticles. In the minimal version of supersymmetry and some of its variants, R-parity 
is assumed to be conserved ad hoc, which prevents single creation or annihilation of superparti- 
cles. However, models with broken R-parity can be constructed naturally, and such models have 
a number of interesting phenomenological consequences. The crucial point is that unlike most ex- 
tensions of the SM, RPV contributes to B decay amplitudes at the tree level. Moreover, the current 
bounds [14] on sparticle masses and couplings leave open the possiblity that such contributions 
can indeed be comparable to or even larger than the SM amplitude. It may be noted that the 
presence of two interfering amplitudes of comparable magnitude is essential for a large deviation 
of CP asymmetries from the SM prediction. 

It is well known that in order to avoid rapid proton decay one cannot have both lepton number 
and baryon number violating RPV model, and we shall work with a lepton number violating one. 
This leads to slepton/sneutrino mediated B decays. Since the current lower bound on the slepton 
mass is weaker than that on squark mass, larger effects within the reach of current round of 
experiments are more probable in this scenario. We start with the superpotential 

W V = KjkLiQjDl (17) 

where k = 1,2,3 are quark and lepton generation indices; L and Q are the SJ7(2)-doublet 
lepton and quark superfields and D c is the £?7(2)-singlet down-type quark superfield respectively. 
This leads to a four-Fermi hamiltonian relevant for B decays [15]: 

H 9 = l d fkn(dnl ti PLdj)s(d k ^P R b) s 

+^df kn (d n j lx P L b)8{d k j lx P R d j )8 

+\uf nk (u nl ' I P L u j )8(d k ^P R b) 8 + H.c. (18) 

where 

A' A'* V \'* A' A'* 

jR _ \ " A »jfc A in3 ,L _ \ ~* A i3k A inj R _ \ - A ijk A in3 ,^ Q , 

a ^~2^ 2m 2 > Sfcn-2y 2ml ' ]nk ~^ 2ml ' 

l u Li I u Li I e Li 
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and the subscript 8 indicates that the currents are in color SU(3) octet-octet combination. 

Following the standard practice we shall assume that the RPV couplings are hierarchical 
i.e., only one combination of the couplings is numerically significant. Let us note that both the 
transitions B — ► rf )_KT(*) and B — ► <f>K are controlled by the quark-level transitions b — ► sss. 
Thus, let us assume, to start with, only c?22 2 an d ^222 to be nonzero, as has been done in [18]. Of 
course, r\K modes can be fed by b — > uus and b — > da's transitions. Since they affect other decay 
modes like 5 — > 7rif where there is no apparent discrepancy with SM expectations, we assume 
those operators to be vanishing. 

Next let us discard d^ too. The reason for this is that d§2 2 an d ^222 are related by SU(2) 
symmetry, and are the same if we neglect the electroweak D-term that causes the sneutrino-slepton 
mass splitting (on the other hand, 0(222 an d ^222 are completely unrelated, and unless there is some 
underlying texture in the RPV couplings, there is no reason why they should be equal). However, 
presence of generates b — > ccs transition, which in turn affects the modes like B — ► J/ipKs 
which is used as a standard to extract sin(2/3). Since the values of sin(2/3) extracted from different 
charmonium modes, as well as from the J/tpir mode with a different quark-level process, are 
almost the same [5], it is a safe assumption, also compatible with the principle of Occam's razor, 
to have no RPV contribution to that channel. Thus, the value of (5 extracted from B — > J/ipKs can 
be taken to be the SM value for that angle. The product coupling d\ 2 2 does not contribute to the 
B° — B° box, so that there is no scope to have an extra box amplitude, in contrast to the situation, 
e.g., in [10] (but d 22 2 contributes to the B s box; this is discussed later). The QCD corrections 
are easy to implement: the short-distance QCD corrections enhance the (S — P) x (S + P) RPV 
operator by approximately a factor of 2 while running from the slepton mass scale (assumed to be 
at 100 GeV) to m h [27]. 

The RPV amplitude for B — ► r/'K is given by 



rfK 



^ "222 



R 



K 



AU 
A T)'K 



1 A* ' 



(20) 



and that for B — > i]K* is 



M„ 
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r,K* — 



-d L 

^"222 



(21) 



where 



Ri = 



m s (m b - m s ) 
,2 



R2 = - 



m 



v 



m s (m b + m s ) 



(22) 



and 



aU,S 

rfK 



fU,S rpt 



AU,S 



For B 



(m v ,)( mi 
2f^m K ,A^ K \m 2 v )(e K *.p v ). 

(f)K (both neutral and charged channels), the RPV amplitude is 

a _ 1 



M 



4>K m 



d222^(t>K 



(23) 



(24) 
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with 



AcpK = 2f (t> m (t) F^ K \m\){e^,.pK) 



(25) 



All these amplitudes are calculated at the slepton mass scale, and, as stated earlier, should be 
multiplied roughly by a factor of 2 when we compute their effects at m^. 

The product coupling ^22 can m general be complex, which we write as 



In our analysis we vary this phase over the range to n, and include the effects of tt < (f)jfc < 2ir 
by allowing jc^l to take both positive and negative values. 

For generic B — > rjK modes, even the SM has two factorizable amplitudes, tree and penguin. 
Following the color transparency argument which predicts the strong phase difference ASsm be- 
tween them to be small [9, 28], we take ASsm = 0. We vary the strong phase difference AS 
between the SM and the RPV amplitudes over a range of to 2ir, but expect to find solutions 
allowing A<5 to be near or 2ir, which is theoretically pleasing. Note that if all strong phase 
differences are exactly zero, there should not be any direct CP violation; indeed, with color trans- 
parency expectations, we should get a small Afip, which anyway is perfectly allowed by data. For 
simplicity we take AS to be the same for all channels to be considered. 

4 The Analysis 

Our input parameters are specified in Sections 2 and 3. We scan the CKM element V u b, which 
has an almost 30% uncertainty, over its entire range. We also vary 1/N C from O.f to f. The 
weak phase <j}£ associated with is scanned over to 7r and the strong phase difference AS 
between SM (tree or penguin) and RPV over to 2tt (we, however, present our results for the 
range — 7r/6 < AS < it/6, motivated by the color-transparency argument). The CKM angle 7 is 
varied between and ir — [3 where (3 = 0.5 axcsm(Acp(B — > J/ifjKs))- 

The following constraints were applied: (i) BR for the modes 7]'K + , r/'K°, r/K* + , 7]K*° and 
4>Ks, (ii) the direct CP asymmetry for — > rj'K^ from BaBar, (iii) the direct CP asymmetry 
for B — > 4>Ks from Belle, and (iv) the average value of BaBar and Belle for sin(2/3) extracted from 
B — ► <f)Ks- Constraints (ii) and (iii) are applied with the rationale that they have larger error 
bars and we wish to check whether the data from the other experiment with smaller errors can be 
accomodated. The direct and mixing-induced CP asymmetries for B — > rj ' K$ are not imposed as 
constraints but one can easily check from the figures that most of the allowed region is perfectly 
compatible with the data. 

Our results are shown in figures 1-4. Let us note the salient features of the analysis. 

• It is known that the nonfactorizable effects in these decays can be substantial. To account 
for that, we have taken N c as a free parameter, and not stuck to its QCD value of 3. 
However, the Wilson coefficients are evaluated with N c = 3. It appears that there is a 
significant nonfactorizable contribution in B — > <f>K, since we have obtained the fit only for 
0.15 < 1/N C < 0.25. Note that r\K channels are N c stable, so the constraint only comes 
from the <j)K mode. This is in conformity with the analysis in [21]. 

• There are two possible bands of solutions, as can be seen from Fig. 1. The left-hand side 
band, with more points and more width, is for negative values of S^ 22 and <j>£ in the first 



^222 = \d222\ ex V{i<t>IjL)- 



(26) 
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Figure 1: The allowed parameter space for A- 32 A- 22 and <f>j[. For more details, see text. 




" 0- -25 -20 -15 -10 -5 5 10 15 20 25 30 
RPV product coupling x 10 4 

Figure 2: One can indeed have solutions with small difference in strong phases of the SM and the RPV 
amplitudes. 
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Figure 3: The range of sin(2/3) as extracted from the B — * <pKs decay. 
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Figure 4: The mixing-induced CP asymmetry for B — > r(Ks- Note that upper half of Band II is disallowed 
from the Belle data. 
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quadrant (this we will call Band I). There is a second narrow band (Band II) for positive 
values of d 22 2 and (f>$. in the second quadrant. This shape is essentially controlled by the 
BRs and CP asymetries of different r\K channels. 

We found a much stronger constraint on the product coupling [A^A'^h 

|A: 32 A'*22l < 2.3 x 10- 3 . (27) 

If we take all these experimental data seriously, it is possible to get even an upper bound: 

lA^A^I > 1.3 x 10- 3 . (28) 

The present bound, as quoted in [29] for the third slepton generation, is only 0.23, and that 
too assuming squarks at 100 GeV. If one has 300 GeV squarks, the bound gets weaker by 
a factor of 9. However, one gets a comparatively better constraint from B s — B s mixing, 
since nonzero values of X' i32 and \' i22 can generate a second amplitude for the box, with two 
sneutrinos and two right-handed strange quarks flowing inside the box. Taking Am# s to 
be completely saturated by the RPV contribution, and using the experimental lower limit, 
one gets a bound on the product coupling which is approximately 1.5 x 10 -2 . The main 
loophole in the analysis is the fact that only an experimental upper bound can generate an 
upper bound on some unknown parameter; moreover, it is questionable to neglect the SM 
contribution completely. Even then we have an improvement by an order of magnitude. 

Figure 2 shows that AS should be positive for Band I and negative for Band II. It can 
be checked easily that this ensures a positive direct CP asymmetry in all channels, and a 
negative mixing- induced CP asymmetry in B — > <pK$- Figure 3 shows the range of sin(2/3) 
in this model, which is between the upper limit of 0.02 and —0.6. The central value of Belle, 
however, cannot be reproduced. 

Figure 4 shows the mixing-induced CP asymmetry for B — ► rj'Ks- Note that imposition of 
the Belle data means a significant portion (the region in the upper right-hand corner) of the 
already weak Band II is ruled out, while Band I is completely allowed. Still, one cannot rule 
out Band II completely. 

The CP asymmetry in B ± — > 7/ K ± is found to lie between 0.115 and 0.22 (the upper limit) 
for Band I, and between 0.14 and 0.22 for Band II. The direct CP asymmetry in the neutral 
channel lies between 0.13 to 0.24. This is definitely compatible with the Belle data, but 
again the central value is far away. One must wait for the error bars to come down. 

The direct CP asymmetry for B — ► 4>Ks lies between 0.13 and 0.25 for Band I and between 
0.13 and 0.22 for Band II. This is in perfect harmony with the BaBar data too. 

Let us try to understand why we get a nonzero Aq P . This is due to the constraint put by 
Aqp (B — ► 4>K). To see this qualitatively, let us assume that there is only one SM amplitude 
and both weak and strong phase differences vanish so that there is no direct CP asymmetry. 
It is easy to see that in that case there is no change in the prediction for sin(2/3), which is 
given only by the phase in the B° — B° box. 

The angle 7 can take any value upto ir — (3 for Band I. On the other hand, it can be only in 
the second quadrant for Band II. The origin of such a pattern can easily be traced back to the 
expressions of BRs and CP asymmetries of the rjK modes. Ref. [22] estimates 40° < 7 < 78° 
at 95% CL from different B — ► tttt and B — > Kir modes. These modes are not affected 
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by c?222! so one concludes that if there are no other nonzero RPV couplings, Band II is 
completely ruled out. Thus, if the value of sin(2/3) extracted from B — > rj'Ks converges 
towards that found from B — > J/ipKg, this solution will be in trouble. However, the error 
bars are too large to draw any definite conclusion right now. 

• As a check, we redo the analysis switching the B — > <pK BR constraint off. There is no sub- 
stantial change in the result, except that the lower bound on Acp{B^ — > rj'K^) marginally 
decreases to 0.11. Thus we consider our result to be fairly robust with respect to hadronic 
uncertainties. However, there is no upper limit on 1/N C anymore. 

We have also explicitly checked that BRs of the so far unobserved r]K channels remain below 
their experimental upper bounds. 

5 Summary and Conclusion 

We found that a minimal set of RPV couplings, compatible with all present data, can explain the 
BR anomalies for B — > rj'K and B — ► r]K*, and the unusual value of sin(2/3) from B — ► ^-ffg. 
This is mainly due to the fact that RPV contributes to B decays at tree-level. A nice feature is 
that one can have points where the strong phase difference between RPV and SM is small, which 
is what one expects from the color transparency argument. 

With more data pouring in, one can significantly shrink the allowed parameter space for 
RPV. However, even at present we get sufficiently strong bound on the relevant product coupling 
— better by two orders of magnitude at least. 

What other effects are mediated by the product coupling A- 32 A- 2 2? This generates a top decay 
channel t — ► ess, whose strength is, unfortunately, only a few per cent of the corresponding SM 
channel. Thus, the signal is essentially unobservable. A better signal may come from the strange 
squark mediated semileptonic top decays t — ► cl + i~ (i = e, //, r). This also generates b — ► svv 
or B s — ► vv decays. The third effect is a new box amplitude in B s — B s mixing. Here, also, the 
amplitude is only at a few per cent level compared to the SM amplitude, so we do not expect any 
significant CP asymmetry in channels like B s —> J/ipcj) even if <j>£ is large. 

We have not discussed the product coupling ci^ which can also explain the anomalies that are 
studied in this paper. As we mentioned before, that coupling also generates the SU(2) conjugate 
transition b — > ccs. This may jeopardize the predictions from the B — ► J/ipKs channel, for 
example. However, a nonzero phase in should show up in the B s box, i.e., one may get a 
significant CP-asymmetry in the channel B s — > J/ip(f>, contrary to the SM expectation. The reason 
is simple: RPV contributes to both the box amplitude and to the decay b — ► ccs. The rare decays 
b — > s£ + £~ or B s — ► £ + l~ also receive a tree-level contribution from an d may be pushed up 
to the observable level. 

At present, one has to wait for the errors to come down. This may rule out part or all of one 
or both bands (existence of Band II is already under threat from the fit of 7, as we have shown). 
However, if the bands still remain allowed, one should look for any unexpected CP asymmetry 
signal in B s decays. Only such correlated studies can unravel the exact nature of new physics. 
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